Reversing quantum dynamics with near-optimal quantum and classical fidelity 
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We consider the problem of reversing quantum dynamics, with the goal of preserving an ini- 
tial state's quantum entanglement or classical correlation with a reference system. We exhibit an 
approximate reversal operation, adapted to the initial density operator and the "noise" dynamics 
to be reversed. We show that its error in preserving either quantum or classical information is no 
more than twice that of the optimal reversal operation. Applications to quantum algorithms and 
information transmission are discussed. 



PACS numbers: 03.67-a, 03.67.Hk, 5.40. Ca, 89.70.+C 

Introduction. Counteracting the effects on quantum sys- 
tems of noise generated by interaction with an envi- 
ronment is a central problem for the emerging field of 
quantum information processing. It's solution can be 
expected to have applications in quantum computation, 
precision measurement and information transmission. In 
this Letter we exhibit a reversal operation which takes ac- 
count of both the noise and the initial density operator, 
to achieve near-optimal preservation of the initial density 
operator's quantum entanglement or classical correlation 
with a reference system. 

We model quantum noise in a system Q by the most 
general dynamics that can arise via a unitary interaction 
{/Q-^ with an environment E in initial state |0^). This 
^ is a trace-preserving completely positive map A : p ^ 
O-i. where A, := {i^\UQ^\0^). The Ai form a 



decomposition of A, A ^ {^i}- When such a map A acts 
on Q, an entangled state \ipo^) ■= X^i v^l^^) 1^*^) Q 
with a reference system R evolves as: 



i 

l*/> ■■=J2VF^\^"')A,\^Q)\.f 



The entanglement fidelity Ff,{p,A) is defined as 
||Po|*/)||2, where Pq := \i^o^){i^o'^\ (S> l"" §. Thus Fe 
is the squared norm of the projection of the final state in 
(|l|) onto the subspace associated with the initial entan- 
gled state IV'o^'^)- It depends only on p :— 
and is given by J^i l^i' ^iPp- When p — \ip){tp\, it is equal 
to the input-output fidelity {■ip\A{\^p){ij\)\ip) . For an en- 
semble E = {pi, Pi} (where state pi occurs with proba- 
bility Pi), we define the average entanglement fidelity by 
Fe{E,A) := ^^^PiFeipijA). A special case is 

F,i{p,A) ■.= J2pMA{m)\t) ^F,i{p,,m},A) , (2) 

i 

where the |i) form an eigenbasis of p and p = 

This is the classical fidelity for the classical information 



(quantified by S{p)) of the ensemble of orthogonal eigen- 
states of the input density operator p. Another special 
case is an ensemble consisting of a single density operator 
p. In this case it is just Fe{p, A) and may be viewed as the 
fidelity for transmission of the amount S{p) of quantum 
information. 

The average entanglement fidelity can equivalently be 
defined as the norm squared of the projection of the over- 
all final state onto the subspace in which entangled states 
I'lpi^'^) representing the initial ensemble are correctly cor- 
related with orthogonal states of an additional reference 
system S: 



F, 



\Pr 



rRS 



RQS\ 



(3) 



where Pc := K^'>(*^'l ® IV^/*'^)(V'f'^l- The initial state 
of RQS may be |i^o^^^) := E^VP'^\l^)\^!^^^)' ^i^h the 
ensemble of entangled states in RQ produced by entan- 
glement with S. Alternatively, the initial state of RQS 
may be mixed, with perfect classical correlation rather 
than entanglement, between a basis of S and the differ- 
(1) ent entangled states of RQ. In this case \'ipo'^^){'4'o'^^ \ is 

(||). The average 



given by E^P^\i ){^ \ m 
entanglement measure is insensitive to whether entan- 
glement, or merely classical correlation, exists between S 
and RQ. 

For example, consider the special case of Fd. Here 
the reference system R plays no role as the pi are pure. 
After suppressing R, Pc = Z^i ^ con- 

tains a record of the classical information sent. S and 
Q may be supposed to be either entangled or classically 



correlated. 



with Pq'^ = 



or 



— SiPir )(* I ® for orthonormal system 

and reference bases where the system basis is the 
eigenbasis of p. In either case, computing the probability 
tr {PcP^'^) that the final system-reference state falls into 
the subspace in which system and reference exhibit per- 
fect classical correlation in the desired bases, gives the 
classical fidelity Fci{p,A). 
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The reversal operation. We motivate the definition of 
the near-optimal reversal operation 7?.^.p by considering 
operations A that are perfectly reversible on a "code" 
subspace C. Let Pc be the projector onto C. Per- 
fectly reversible operations have a decomposition Ai for 
which AiPc = ^ypiWi, where the Wi are partial isome- 
trics from C into orthogonal subspaces, which means that 
W^Wj = S^jPc m. Without loss of generality, as- 
smxie that the ranges of the Wi span the state space. 
The reversal operation has a decomposition consisting of 
the operators = PcA] / y/pi. This resembles the ad- 
joint (defined using the Hilbert-Schmidt inner product 
{A, B) := tv A^ B on operators) of the restriction Ac 
of A to C, which is given by y^J. ~ {PcA\} = {^iW}}. 
The ^Jpi need to be removed to ensure that the reversal 
operation is trace preserving. The general definition of 
'R-A,p is also based on the adjoint, suitably corrected to 
ensure that it is trace preserving (and continuous in the 
density operator): 

nj,,,^{p''^A\A{p)-''^}. (4) 

This agrees with the reversal operation for codes given 
earlier, using the uniform input state p — Pc/i^^ Pc- 

The notation TZa.p is justified by: 
Lemma 1. The definition of the reversal operation TZa.p 
is independent of the decomposition {Ai\ of A. 
Proof. Let A ^ {Bi] be another decomposition. By 
adding null operators to one of the two decompositions, 
we can ensure that both have the same number of op- 
erators. Note that adding null operators to {Ai} does 
not change the action of TZa,p- Then there exist Uij such 
that Bi = UijAj, where the matrix u defined by Uij 
is unitary. The decomposition of 'R'A,p given in trans- 
forms via the coefficients of into a decomposition given 
in terms of the Bj. As is also unitary, the result is 
another decomposition of the same operation. ■ 

The operation TZa,p is near-optimal in the sense given 
by the following theorem. 

Theorem 2. Let E = {pi , pi} be an ensemble 

of commuting density matrices, and let p J^iPiPi- 
Then for any trace-preserving completely positive map 
7^, F,{E,TZa,pA) >F,{E,TZA)\ 

As a corollary, if Fe{E,TlA) = 1 — ?/, then 
Fe{E,TZA.,pA) > (1 - ry)^ > 1 - 2??. That is, 7^^,p's 
error is never greater than twice that of the best reversal 
operation. 

Proof. Without loss of generality, assume that TZ ^ 
{i?i}'s domain is the algebra of operators on supp(.A(/9)) 
and its range is supp(/9). Allowing more general reversal 
operations cannot increase entanglement fidelity. Then 
there exist operators Bi such that 

i?. = p'/'BlA{pr'/' , (5) 

namely those defined by i?| := p~^^^RiA{py^^ . (Gener- 
alized inverses are to be understood here.) Let B ~ {Bi}. 



We have 

Fe(i^,7^^) =^PZ^|trpl/2st_4(^)-l/2^^.^^|2 (g) 
I ij 

Define Xlj := tr pi/2^|^(p)-i/2^]p/. By proper 
(Z-dependent) choice of operator decompositions B ~ 
{B-} and A ~ {Al} (corresponding to singular value 
decompositions — cf. §7.3 — of the matrices X'), we 
may obtain the same expression, but with the inner sum 
having just one index. Then applying the cyclicity of 
the trace, the fact that [p, pi] ~ 0, and defining Xu := 
p\l'A{pr^l^A\p^l^p\l\ Yu ■.^p\l'A{p)-^I^B\p^l^p\l\ 
gives: 

F^{E^TIA) ^Y^ViY. P^'^b'^A{p)-^I^A\pA^ 

I i 
il il 

<{Y,\trXlXu\'Y.\^rY^,YM'Y/^ 

il i'l' 

<(^|trXtx.nV2<(^|trXtx,,f)V2 

il ijl 
I ij 

= Fe{E,nA,pA)^^^ . (7) 

Here the first two inequalities are Schwarz inequalities, 
the third uses the fact that 

^|trpi/2st_4(^)-i/2^^^^|2 < F^^p^TZB) < 1 , (8) 

i 

the fourth just adds positive terms inside the square root. 
The last identity depends on Lemma |l|. Since B is not 
necessarily trace-preserving, (|^) is not immediate, but 
from (|^) and the trace-preserving condition on TZ, 

^ RlR, = A{p)-'/'B,pBlA{p)-'^' = 1. (9) 

i i 

It follows that B{p) = A{p), a normalized density oper- 
ator. Let iV'o''^) be a purification of p. Then the states 
(I«®i3«)|7/'^^)(V'^^| and (X«®7^'^e'3)|^«Q)(^«Q| are 
also normalized density matrices, whence Eg (p, TZB) = 
tr Poil'^ (g> 7^QS«)|V'?'^)(^/'^'5| < L ■ 

Important special cases of Theorem ^ are noted in: 
Corollary 3. For any trace-preserving completely posi- 
tWeU^Fdip^A) < y^Eciip, nA,pA) and Fe{p,nA) < 
^Ee{p,nA,pA) . 

When the members of the input ensemble pi do not 
commute, we do not know whether TZa,p, for p := 
^iPiPij is still near-optimal. 

Relationship to the "pretty good measurement". The 
above analysis of the fidelity of reversal makes it clear 
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that Ti-A,!) provides a method for distinguishing, with 
close to optimal average error, density matrices from 
the ensemble {pj,/5j}, where pj :— A{\i){j\), and p = 
Pj 1 ■ This provides a near-optimal method for dis- 
tinguishing density matrices in an arbitrary ensemble, for 
any ensemble {Pj,Pj} may be constructed by an opera- 
tion 



(10) 



where pj = 



are the spectral decomposi- 
tions of the density matrices to be distinguished. The 
operation A may be thought of as measuring in the or- 
thogonal basis and then producing the corresponding 
/5j, for example by randomly applying, with probabilities 
Aij, unitary rotations taking |j) into With A as 

defined above. 



-1/2. 
•ut J 



(11) 



The "pretty good measurement" (PGM) was introduced 
by Holevo [g (the term "pretty good measurement" is 
from 1^) for the case of linearly independent pure states, 
in which case the PGM is a measurement of orthogo- 
nal projectors, and as Holevo showed, the optimal such 
measurement. For an ensemble of unnormalized density 
matrices pj := PjPj, where Pottt ■= J2j Pj ^ normalized 
density operator, the PGM is defined by the set of oper- 
ators consisting of the 



Y -1/2 -1/2 

^j Pout PjPout 



(12) 



For pure pj cx these are just the operators corre- 

sponding to the "/9-distorted" states Pout^lj)- Note 
that for a doubly indexed ensemble of unnormalized 
states pij, J2i-^ij — ^j I where the Xj are the PGM 
for the Pj := pij. The operation ( pi] ) may be viewed, 
via the given representation, as performing the PGM 



for the ensemble of unnormalized states ^ypj\/\j\vij), 
and returning |j) when the measurement result ij is ob- 
tained. Indeed 
therefore J2i ^h^ij 



for that ensemble RljRij 



Xij, and 



Xj. Thus the operation may also 
be viewed as doing the PGM for the pj, and return- 
ing \j){j\ when the measurement result is j. However, 
a given ensemble pj may in general arise from orthogo- 
nal states \j) by actions of channels different from the 
"classicizing" one (|lO|), which completely decoheres the 
orthogonal states \j) before producing pj. For example, 
if the Pj are orthogonal and pure they may be produced 
either by measurement in the basis |j) followed by an ap- 
propriate unitary operator U, or by applying U without 
prior measurement. In the first case quantum coherence 
is completely destroyed, while in the second case it is 
perfectly preserved. When the channel producing the pj 
is not of the form dlQ), the reversal operation TZ^.p will 
be different from (pT[). Although (|l^) still gives near- 
optimal classical fidelity, it will not necessarily give good 



entanglement fidelity, since in some sense it decoheres 
the states pj. 'R-a,p: however, will have near-optimal en- 
tanglement fidelity while retaining near-optimal classical 
fidelity. 'R-A.p thus takes advantage of whatever coher- 
ence remains between the pj ; it avoids decohering the pj 
if the channel has not decohered them already. 
A hound on the classical fidelity of reversal. To bound 
a fidelity of reversal it is sufficient to bound the fidelity 
for the near optimal reversal operation and apply The- 
orem |[ Here we have a look at such bounds for classi- 
cal fidelities of reversal for A of the form (10). In this 



case, the classical fidelities are probabilities of success for 
measurements that attempt to infer which of the \ of 
the input state p = ^jPj\j){j\ was actually transmitted. 
The expression for the PGM gives the following bound 
on the optimal probability of success Fd (with the defi- 
nitions of the previous section): 

< F,i{p,nA,pA) 

E, , -1/2 -1/2 N 
tr KPout PoPout Pj) 



1 - XI ^PouC^PiPouPpj), 
i-.j-.i^'j 



(13) 



where we used the identity E»j iPoul'^ PiPouP Pj) 
1. Thus the probability of error E^i is bounded by 

tr {PouI'^PiPouI^P])^ which is a multiple of the 
sum of the Hilbert-Schmidt inner products of the differ- 
ent Pout^ PjPout^ ■ When Pout is uniform (proportional to 
a projection), this sum can be easy to estimate. An often 
used measure of overlap between density matrices is the 
Bures-Uhlmann fidelity . This measure depends only 
on the pair of density matrices, not on an overall aver- 
age, and is defined by Fbu{<^i,'^2) ■= tr 



1/2 1/2 



cr2Cri 



The expression for the optimal reversal given in (11) can 
be used to derive a bound on the probability of error in 
terms of the Bures-Uhlmann fidelities. 
Theorem 4- 



F^i > 1 - X! VP^^BuiPiJj) 



i,j-i¥'j 



(14) 
(15) 



in a row. 



Let Aj be the matrix whose i'th column is 
I , and A the matrix obtained by placing the 
Then AA^ — pout- Let Rj be the matrix 



Proof. 

A 

1/2 

whose j'th row is ^JWjy ^ij{''^ij\Pout ^^'^ R the matrix 
obtained by placing the Rj in a column. R is simply an 
explicit matrix form of TZa,p- Since R = yl^(ylAt)^i/^, 
i? is a matrix with the property that RA is positive semi- 
definite. (In fact, this gives an alternative approach to 
defining TZa,p)- The matrix RA has a natural block struc- 
ture that mirrors that used to define R and A. It is read- 
ily verified that Fd for TZa,p is given by the sum of the 
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squared Frobenius norms |(i?^)jj|2 := tr j^- 
of the diagonal blocks. Since \RA\l is one, it suffices to 
estimate the sum of the squared Frobenius norms of the 
off-diagonal blocks to bound the optimal Fd- To do so, 
observe that (RA)"^ = A. The block at block posi- 
tion l,k in B :— A"^ A is given by Bik — Ai''Ak. Since 
AiAi'' = pi, 



Bik^'Bik = Ak^'piAk 



jj 1/2 l/^jji 

Upk PiPk 



(16) 



(by polar decomposition of A^) for some unitary oper- 
ator U. Consequently, the ii-norm of Bik, defined by 

\Bik\i ■■= tr a/ Bik^Bik is given by y/piPkFBu{Ph Pk)- It 
therefore suffices to relate the Frobenius norms of the off- 
diagonal blocks of a positive semi-definite matrix to the 
Li norms of the off-diagonal blocks of its square. 

a 

Lemma 5. Let M = \ , | be positive semi-definite. 



with a, 6, c matrices. Write Af 



y z 



with the 



same block structure. Then < \b\i. 

Proof. Without loss of generality, assume that y is 

non-negative diagonal. Otherwise, with a block-diagonal 

unitary [/ = ^ ^ ^ ^ with u and v chosen to implement 

the singular value decomposition of y, we may transform 
M and M^/^ so that y is a diagonal matrix with non- 
negative diagonal entries (for rectangular y, the upper 
or left-hand square portion is diagonalized). This does 
not affect the norms, since U transforms blocks indepen- 
dently and the Li and Frobenius norms are both unitar- 
ily invariant. Let yi be the diagonal entries of y. Note 
that b ^ yx + zy and > tr b (cf. g, p. 432). Now 
tr {yx + zy) = J2zyii^ii + ^a)- By the positivity of M^/^^ 
Vi < XiiZii, so yi is less than at least one of Xu, zu. Thus 
lyli = Ei yf < tr 6 < as desired. ■ 
To apply Lemma ^, consider first the block decom- 
position of RA and B determined by Bn. The squared 
Frobenius norm of the first block row and column exclud- 
ing [RA)ii is bounded by the sum of the Li norms of the 
corresponding block row and column in Bn. By subad- 
ditivity of the norm, this is at most X)i>i(l^ii|i + l^iiD- 
After a suitable permutation, the same argument applies 
to the the row and column determined by Bn, for each 
i. The proof of the theorem then follows by summing 
over the resulting inequalities and noting that each off- 
diagonal block occurs twice on both sides. ■ 
Applications. Due to its near optimality, the reversal op- 
eration Ti-A,p can be used in any situation where classi- 
cal or quantum information has been corrupted by noise 
with known behavior. 'Ti-A,p has a simple definition, 
but whether it or a good approximation can be imple- 



mented efficiently depends on the details of the situa- 
tion. Whether or not it can be efficiently implemented, 
because its error is at most twice the optimum, it can 
be used as a theoretical tool to obtain upper bounds on 
the achievable fidelities in a given situation. The upper 
bounds can then be compared to the fidelity achieved by 
simpler algorithms. An example of this occurs in the use 
of stabilizer codes for quantum error-correction. When 
the noise model is independent and depolarizing, classi- 
cal coding theory immediately suggests a combinatori- 
ally straightforward error-correction algorithm based on 
maximum likelihood error syndrome decoding. Compar- 
ing this method to 'R-a,p suggests itself as a fruitful path 
of investigation with applications to asymptotic bounds 
in quantum coding theory |ic| ]. 

Another application is to query complexity for quan- 
tum oracles. Here we are given a quantum black box im- 
plementing an unknown quantum operation from some 
set. A simple method for attempting to determine which 
operation we are given is to apply it to copies of some in- 
put state and attempt to distinguish the output state. A 
bound on the probability of success can then be obtained 
by using bounds such as the one of Theorem |[ This was 
how the fact that the hidden subgroup problem has low 
query complexity was first realized . 
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